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Abstract. We describe an optical fiber based interferometer to measure velocity profiles in sheared com- 
plex fluids using Dynamic Light Scattering (DLS). After a review of the theoretical problem of DLS under 
shear, a detailed description of the setup is given. We outline the various experimental difficulties induced 
by refraction when using a Couette cell. We also show that homodyne DLS is not well suited to measure 
quantitative velocity profiles in narrow-gap Couette geometries. On the other hand, the heterodyne tech- 
nique allows us to determine the velocity field inside the gap of a Couette cell. All the technical features 
of the setup, namely its spatial resolution (~ 50-100 /im) and its temporal resolution (« 1 s per point, 
~ 1 min per profile) are discussed, as well as the calibration procedure with a Newtonian fluid. As briefly 
shown on oil-in- water emulsions, such a setup permits one to record both velocity profiles and rheological 
data simultaneously. 
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1 Introduction 



The behaviour of complex fluids under flow raises impor- 
tant fundamental problems and involves a very wide spec- 
trum of possible applications. A complex fluid is charac- 
terized by a mesoscopic scale, located somewhere between 
the microscopic size of the molecules and the macroscopic 
size of the sample . For instance, for oil-in- water emul- 
sions, this mesoscopic scale would be the size of an oil 
droplet. The existence of this mesoscopic scale can make 
the flow of a complex fluid considerably difficult to un- 
derstand as compared to that of a simple, viscous fluid. 
Indeed, most of complex fluids show strong viscoelasticity 
and/or important nonlinear rheological behaviours. 

These behaviours arise from coupling effects between 
the fluid microstructure and the flow. Usually, the pres- 
ence of a shear flow tends to order the fluid, for instance, 
by aligning chains in a polymer solution and thus decreas- 
ing the viscosity ("shear thinning"). On the other hand, 
in a class of concentrated systems recently called "soft 
glassy materials" , such as concentrated emulsions or col- 
loidal glasses, yield stress and jamming phenomena are 
frequently encountered [^[Q]. 

Although the characterization of the global flow prop- 
erties remains a first essential step, much interest has 
grown in performing local measurements on complex fluids 
under shear. Indeed, in many cases, inhomogeneous flows 
are observed for which the usual analysis in terms of pure 
shear flow is not applicable and may be misleading. 



This paper is devoted to local velocimctry in sheared 
complex fluids using Dynamic Light Scattering (DLS). 
The need for local measurements is highlighted in the 
next Section through three examples taken from recent 
works in complex fluids. We then review a few existing 
techniques that provide local measurements and explain 
how they may or may not be suited to experiments in 
the concentric cylinder geometry (referred to as "Couette 
geometry" in the following). In Sec. ||, we give the basic 
theoretical principles of DLS under shear. In Sec. f|, we 
present an experimental setup based on the use of single 
mode optical fibers. Finally we show that the heterodyne 
geometry is better suited than the homodyne one when 
the fluid under study is sheared in a Couette cell. 



2 The need for local measurements in the 
rheology of complex fluids 

2.1 Examples of inhomogeneous flows in complex 
fluids 

Complicated rheological behaviours are observed in com- 
plex fluids whose flow curves cr(j) are non-monotonic (a 
stands for the shear stress and 7 for the shear rate) . In such 
cases, above a critical shear rate, the stress may decrease 
resulting in a mechanical instability of the flow ||. For 
controlled-shear-rate experiments, this type of instability 
is expected to show up on the flow curve as a stress plateau 
for 71 < 7 < 72 separating two regions where the stress is 
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an increasing function of the shear rate. Such a signature is 
indeed observed for various complex fluids such as worm- 
like micelle solutions Q or some lyotropic lamellar phases 
0. On the stress plateau, the fluid is believed to phase- 
separate into two regions: one with high viscosity and low 
shear rate flowing at 71 and one with low viscosity and 
high shear rate 72. This phenomenon is known as shear 
banding. In micellar solutions, each "band" was shown to 
correspond to different fluid structures, isotropic or nc- 
matic, using flow birefringence || or NMR spectroscopy 
[^1 . The correspondence between the branch flowing at 72 
and the nematic state is still under debate |J . 

Thus, even though the geometry and the nature of 
the deformation imposed on a complex fluid remain fairly 
simple (constant, homogeneous force or velocity at the 
wall of a Couette cell for instance), the response of the 
fluid may be highly nonlinear and sometimes spatially in- 
homogeneous. In some extreme cases, one not only ob- 
serves bands of different structures and flow properties, 
but also large temporal fluctuations of the global sample 
viscosity M,OJl2|. This last behaviour, recently coined 
"rheochaos" uW, is observed in the vicinity of the shear- 
induced transitions described above. The occurence of such 
dynamical behaviours shows that inhomogeneous flows of 
complex fluids can exist both spatially and temporally un- 
der simple shear, even in the absence of strong inertial or 
elastic effects [|l4]|n|. 

A simpler but recurrent problem in rheology is caused 
by wall slip. Unlike simple fluids, mixtures under shear 
may involve thin lubricating layers of low-viscosity fluid 
in the vicinity of the cell walls. These lubricating films 
absorb most of the viscous dissipation so that the fluid in 
the bulk is submitted to a shear rate smaller than that im- 
posed globally. In emulsions, it is generally believed that 
apparent slip is due to "wall depletion" , the fluid close 
to the walls containing less oil droplets than in the bulk 
PI . In classical rheology, wall slip is usually investigated 
by varying the cell characteristics and/or the wall rugos- 
ity Q. In any case, it remains a problem of considerable 
importance for industrial applications and shows yet an- 
other occurence of inhomogeneous flows in simple shear 
geometries. 

2.2 Experimental constraints in the Couette geometry 

Through the three examples of shear-banding, "rhcochaos,' 
and wall slip, it clearly appears that local measurements 
are crucial to understand the physics underlying the re- 
sponse of many complex fluids to deformation and flow. 
However, most of the studies in this area still heavily rely 
on global data recorded by standard rheometers imposing 
a torque (a velocity resp.) on the axis of a moving cylindri- 
cal or conical part and measuring its velocity (the torque 
resp.) once the fluid is set into motion. 

For more than two decades, great experimental effort 
has been devoted to dcvclopping new methods for mea- 
suring local quantities of rhcological interest such as the 
local stress, the velocity field, or the local velocity gra- 
dient [^|]g| Jl7| . |l8| . Ideally, such a technique should be im- 



plemented on standard rheometers. On these apparatus, 
various geometries are used to induce a pure shear flow: (i) 
Couette cells in which the fluid is confined in the gap e be- 
tween two coaxial cylinders and where the shear rate can 
be considered as homogeneous, provided the ratio of the 
gap to the cylinder radii e/R is not too large, (ii) cone- 
and-plate cells which provide the best approximation of 
pure homogeneous shear flows, and (iii) plate-and-platc 
cells where the shear rate increases linearly with the ra- 
dial position. 

In the present paper, we will restrict ourselves to shear 
flows generated in a Couette cell. This particular geometry 
imposes a few important constraints: 

1. The measurement technique has to be non intrusive. 
Indeed, the very high sensitivity of the material to lo- 
cal deformations precludes the use of a probe, such as 
the hot wires or the hot films classically used in hy- 
drodynamics to record fluid velocities. 

2. The thickness e = 1-3 mm and the optical properties of 
the samples under investigation usually hinders direct 
optical access needed for techniques like Laser Dopplcr 
Velocimetry (LDV) or Fluorescence Correlation Spec- 
troscopy (FCS) : most complex fluids are not transpar- 
ent and strongly scatter light. Also, in general, we are 
not able to directly follow individual droplets in an 
emulsion or fluorescent seeding particles in a lamellar 
phase. Let us point out, however, that direct track- 
ing of the system is possible in the particular case of 
two-dimensional foams [jl9f or when the fluid is trans- 
parent enough to allow Particle Imaging Velocimetry 
(PIV) measurements such as those by Pine et al. in 
Ref. |. 

3. Finally, the curved geometry of a Couette cell induces 
refraction effects that can be relatively strong and make 
optical techniques relying on light scattering tricky to 
implement. 

Moreover, a reliable measurement method should yield 
accurate, reproducible measurements with good temporal 
and spatial resolutions. Here, we shall call "good" a spatial 
resolution which provides, say, 20 measurement points ac- 
cross the gap of our Couette cell, which is typically e = 1- 
3 mm wide. Any technique that is able to distinguish vol- 
umes of fluid separated by 100 (im and measure their re- 
spective velocities (or any other local quantity) will thus 
have a good resolution. 

2.3 NMR and DLS, two local measurements 
techniques 

Among the few quantitative local techniques that have 
proved successful for complex fluids under flow are (i) Nu- 
clear Magnetic Resonance (NMR) and (ii) Dynamic Light 
Scattering (DLS). 

2.3.1 Nuclear Magnetic Resonance 

Callaghan et al. have used NMR to image the velocity 
field of sheared wormlikc micelles in both cone-and-platc 
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and Couette (e = 1 mm) geometries Results in the 

isotropic-nematic coexistence regime show a strong corre- 
lation between shear bands and bands of different struc- 
tural order. NMR can resolve details down to 20 (im J|,|)| 
and appears as a very powerful tool for local measure- 
ments in complex fluids. However, the main drawback of 
this technique is its poor temporal resolution: in Rcf. Q 
the authors report on a broad velocity distribution due to 
fluctuations of the flow field but are not able to resolve 
these fluctuations in time. Indeed, the accumulation time 
needed for one profile is about an hour. 

Moreover, this technique requires the costly produc- 
tion and use of strong magnetic fields. Also, the flow cells 
have to be made nonmagnetic. This leads to even more 
important constraints on the design of the experiment. 



2.3.2 Dynamic Light Scattering 



Although subject to restrictions (2) and (3) listed above, 
DLS is a low-cost technique easier to set up, at least in 
the homodync geometry. In the present work, these two 
restrictions can be partially avoided by controlling the op- 
tical index of the complex fluid under study. For instance 
by matching the indices of the aqueous phase and of the 
oil phase in an emulsion, we are able to control the amount 
of scattered light and avoid multiple scattering |^0[ 

DLS is a common tool for probing spatio-temporal cor- 
relations of the fluctuations of the dielectric permittivity 
With a realistic model of these fluctua- 



in a material 21 



tions, DLS allows one to obtain information on the dynam- 
ics of the fluctuations at length scale A, the wavelength 
of incident light. Such fluctuations may arise from differ- 
ent microscopic origins: if scattering particles are subject 
to Brownian motion then DLS yields a measurement of 
the diffusion coefficient of these particles 21 ; in lyotropic 



lamellar systems, fluctuations arise from undulating mem- 
branes and DLS may be a useful tool to characterize the 
elasticity of the surfactant bilaycrs . 

In the simple case of a shear flow, it has been shown 
that homodyne DLS may give information on the local 
shear rate. Such a technique has been used to measure 
shear rate profiles of polymeric fluids in four-roll mill flows 
|l7| , p3| . However, this homodyne technique is not very well 
suited to the Couette geometry as we will see in the follow- 
ing. Ackerson and Clark have shown that heterodyne 
DLS was a powerful technique to measure local velocity 
and applied it to colloidal crystals under Couette flow. 
With the same kind of heterodyne setup, Gollub et al. 
studied velocity profiles in Newtonian fluids at the onset 
of the Taylor-Couette instability |^,^6|. In spite of the 
above results, very few local measurements using hetero- 
dyne DLS have been reported in the rheology literature. 
This is certainly due to the experimental difficulties to set 
up an interferometer around a Couette cell. 



3 Theory of Dynamic Light Scattering in a 
shear flow 

3.1 Homodyne and heterodyne DLS 

The aim of this section is to give the basic theoretical 
principles of DLS under shear which are necessary to re- 
duce the optical information to local shear rate and local 
velocity measurements. In a typical DLS experiment as 
sketched in Fig. |], the sample is illuminated with coher- 
ent polarized laser light. Scattered light is collected at an 
angle 9 by a photomultiplier tube PMT. 




PMT 



Fig. 1. Geometry of homodyne DLS experiments. PMT 
stands for photomultiplier tube. Other symbols are defined in 
the text. 



The optical device which collects light is composed of 
a field diaphragm Df and a lens L. The scattering volume 
V is the intersection of the incident beam with the im- 
age of Df through L. The aperture diaphragm D a allows 
one to reduce the range of collected angles 9 (it controls 
the size of the photodctcction area) . The auto-correlation 
function C(r) =<i(t)i(t + t) > of the scattered intensity 
i(t) is measured by an electronic correlator as a function 
of the time lag r. ej is the direction of the polarization of 
the incident beam and ef is the direction of the analyzer 
between the scattering volume and the PMT. 

3.1.1 Homodyne DLS 

When the dielectric permittivity tensor e(r, t) of the mate- 
rial is slightly perturbated, i.e. when e(r, t) = eT+Se(r, t), 
and Se <C e (I is the unity tensor), the magnitude of the 
electric field scattered at angle 9 is [ElJ: 

E s (t) cx <5ei/(q, t) , 

where q = k/ is the scattering vector, whose mag- 
nitude is given by q — (Ann/X) sin (0/2), and 5ei/(q, t) 
is the « /-component of the Fourier transform of the di- 
electric permittivity tensor of the sample restricted to the 
scattering volume V: 

Se if (q,t)= [ d 3 re 4q r ei ^(r, t)e t . 
Jv 
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In the following, we will focus on the zz-component of 

the tensor Se(r,t), and we write for simplicity: 5e(r,t) = 
8e zz {r, t). 

Measuring the intensity scattered at angle 9 allows 
one to determine the spatial fluctuations of the dielec- 
tric permittivity at length scale <? _1 . In DLS, temporal 
fluctuations are studied through the correlation function 
C(r) =< i(t)i(t + t) >. If <5e(r, i) is a Gaussian random 
variable, one can show that pl|: 



where 



C(r)=A + B\g( 1 \r)\ 



g W(r) =<Se(d,t)5e(c L ,t + T)> 



(1) 



A and B are two constants that depend on the signal-to- 
noise ratio, g^'ir) is the auto-correlation of the Fourier 
transform of the fluctuations of Se(r, t) at length scale q~ l 



and time r. A device measuring 
modyne DLS setup. 



(t)| 2 is called a ho- 



3.1.2 Heterodyne DLS 

In the case where the temporal fluctuations of Se(r, t) are 
non-Gaussian or when g^'ir) contains phase terms, ho- 
modync DLS is not sufficient. To access the real part of 
g^'ir), one has to use DLS in heterodyne mode. This 
method consists in measuring the signal resulting from the 
interference between the scattered electric field E s and a 
local oscillator E LO taken from the incident beam, whose 
magnitude is much greater than the scattered electric field. 
In this case, one can show that the autocorrelation func- 
tion of E s + E LO satisfies pl[: 



C(t) = A' + B'lle (g {1 '(r) 



(2) 



where A' and B' are two constants that depend on the 
signal-to-noise ratio. 

3.2 A simple model for g^'(r) in a shear flow 

3.2.1 Physical considerations 

Let us now turn to the case when the scattering volume is 
submitted to a shear flow, as shown in Fig. The theoret- 
ical problem of DLS in a shear flow has been extensively 
studied The main result of these works is 

that homodyne DLS allows one to measure the local shear 
rate through decorrelation terms of geometrical nature. 

On the other hand, heterodyne DLS yields an estimate 
of the mean velocity in the scattering volume using the 
well-known Doppler effect. Indeed, since scattering parti- 
cles are moving at velocity v, the frequency of the scat- 
tered electric field is Doppler shifted by the quantity qv. It 
is straightforward to show that g^'(r) = F(t) exp(iqv-r), 
where F(t) is some decreasing geometrical factor. Since 
homodyne DLS measures the modulus of g^{r), the lo- 
cal velocity can only be accessed by using the heterodyne 
configuration. 




V 



k; 







a 






\ 



0- 




Fig. 2. Geometry of the scattering volume in a shear flow. V 
stands for the scattering volume and A for the photocathode 
area. 



The aim of this section is to give some simple physi- 
cal ideas that help to understand the correlation functions 
measured with heterodyne DLS. We then present more de- 
tailed equations which will be solved numerically in Sec. 

The geometry of DLS under shear is shown in Fig. 
The incident beam is normal to the velocity (y) and to the 
vorticity (z) directions. In a shear flow, one expects the 
following invariance rule for the correlation of the fluctu- 
ations fl24|: 



r(r, t, r', t')= < 5e(r, t)Se(r', t') > 

= r(r'-r- (v + <y(r.e x )e tf )(t / 



t),t'-t). 



(3) 



Thus, in the presence of a shear flow, the physical prop- 
erties in the scattering volume are no longer invariant un- 
der spatial translations. Hence the spatial coherence of the 
scattered electric field is modified [Q. From Eq. (||), one 
can show that: 



<(5e(q,i)5e(-q',i + r)>^ 0, 
only if the following condition is satisfied: 
q' = q + jr(q.e y )e x . 



(4) 



(5) 



In fact, this relation is strictly true only for an infi- 
nite scattering volume. In practice, V corresponds to the 
intersection of two infinite cylinders. Therefore, the scat- 
tering volume has dimensions of the order of a x a x c (cf. 
Fig. |2|), and the scattering vectors q are only defined up 
to the precision 8q x = ±2tt/c and Sq y = Sq z = ±27r/a. 
Due to the finiteness of the scattering volume, light scat- 
tered at a distance R corresponds to speckle grains and is 
spatially coherent in small coherence areas of typical size 
XR/c. In the experiments, the aperture diaphragm D a col- 
lects a small range of angles 9 and therefore a small range 
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of q's. In this more realistic case, the correlation function 
measured with heterodyne DLS reads: 



9 



(i) 



(r)= J2 <Hq,t)5e(-q',t + T)> 



(0) 



q,q'6-4 



where the sum is calculated over the coherence areas which 
are collected on the PMT. In the presence of shear, due to 
the condition of Eq. (||) , one can see that coherence areas 
are moving in the ^-direction. This leads to dccorrclation 
when ||q' — q|| > Sq x i.e. when r > = 2ir / \cyq y ). 

Another decorrelation time arises from the finiteness 
of the photocathode area. Indeed, speckle grains move out 
of the photodetection area due to the spatial coherence 
condition of Eq. (|^). The separation 5R at a time r be- 
tween two coherence areas which contribute to C(r) is of 
the order of 5R « ||q' — q\\R/kf « •yRr. Thus, for times 
t > T4 = b/(^R), this term leads to the decorrelation of 
C(t) because the two coherence areas are separated by 
more than the characteristic length b of the photocathode 
area. 

Finally there is a third decorrelation time due to the 
transit time r t of the particles through the scattering vol- 
ume V. This transit time is of the order of r t w a/v, where 
v is the mean velocity in the scattering volume. In a lot 
of experimental situations however one has: r t 3> T-y and 
T t 3> ta so that this decorrelation cannot be observed. 
Physical systems may also display intrinsic dynamics such 
as Brownian motion. Such dynamics leads to decorrelation 
but in most systems the decorrelation times associated to 
those dynamics are very large compared to the geometri- 
cal times discussed above and cannot be observed Q . 

From these physical considerations, we can assess that 
<7^(t) takes the following form: 



g^( T ) = e^F(T A ,T^T) 



(7) 



3.2.2 Derivation of F(t^,t^,t) 



We now present a detailed derivation of the heterodyne 
correlation function under shear. In order to estimate C(r) 
in realistic conditions, one must take into account the 
finiteness of the scattering volume V and of the photo- 
cathode area A as sketched in Fig. |[ One can show that 
the magnitude of the electric field scattered at point R = 
R + r' from the scattering volume V is given by plj : 



E s (R,t) = 



k 2 Eoe i{k,Ro-u>t) 

— 1 

47re R 



Se q + 



R„ 



(8) 



with q = kf — ki and assuming that R r and R 3> r' . R 
is the distance between the center O' of the photocathode 
area A and the center O of the scattering volume V. 

In heterodyne DLS, the normalized correlation func- 
tion takes the following form: 



C(r) 



Tie (^JJ^ <E s (R',t)E* 3 (R" 



i + -r)> 



(9) 




Fig. 3. Typical configuration for light scattering in realistic 
conditions. V is the scattering volume, A the photocathode 
area 



Using Eq. (||), one gets: 



C(r) oc Tie J J dr'dr"e ik f< r '- T " 
Ro 



~4-r',t) 8e ( q+ ^r",t + r 
Ro 



(10) 



Including the effect of shear requires a model for the 
dielectric permittivity of the material. A simple approach 
consists in considering isotropic, independent, point-like 
scatterers moving in a shear flow. In such a case, 5e(r, t) = 

J2i=i$( r - r jl£)) an d r i(^ + r ) = r i(*) + vt + jxe y T (see 
Fig. ||). Eq. ( |l0j ) can then be rewritten as : 



C(t) (xlle dr'dr"e lk f {r '- r " } 



dr e it l'( v +T' :re !')' r e i (' c //- R o)[r-(r'-r")+r'-(v+7a;e !/ )T 



(11) 

To proceed further in the discussion, the shape of the scat- 
tering volume and of the photocathode area have to be 
specified, so that Eq. (0) can be solved numerically. 



4 DLS experiments in a Couette flow 

4.1 Description of the experimental device 

As discussed in the Introduction, one of the major chal- 
lenges in the study of complex fluid flows is to measure 
simultaneoulsly the global rheology of a material under 
shear and the velocity profile inside the fluid. In order to 
perform such experiments we have developed the setup 
sketched in Fig. where heterodyne DLS has been in- 
stalled around a classical rheometer (TA Instruments AR 
1000). 

In order to probe the rheological properties of our 
complex fluids, we use a transparent Mooney-Couette cell 
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*-r*SF. 



Fig. 4. Heterodyne DLS setup. La denotes lasers, BS a beam 
splitter, SF spatial filters, and C the device coupling optical 
fibers / 



whose geometrical characteristics are given in the next 
paragraph. The temperature of the sample is controlled 
within ±0.1°C using a water circulation around the cell. 
The rheometer sits on a mechanical table whose displace- 
ments are controlled by a computer. Three mechanical 
actuators allow us to move the rheometer in the x, y, and 
z directions with a precision of 1 /j,ui. 

A He-Ne polarized laser beam La x (40 mW, A = 
632.8 nm) is directed through the cell and is focusscd in- 
side the gap by a lens L (/ = 10 cm). The incident beam is 
polarized along the z direction. A spatial filter SF 1 , com- 
posed of a microscope lens and a single mode optical fiber 
/i (core diameter d = 4 fim, IDIL) collects the electric 
field scattered at an imposed angle 8i by the sample lying 
in the gap. 

In order to perform heterodyne measurements, a lo- 
cal oscillator is obtained from the incident beam by a 
beam splitter BS and a second spatial filter SF 2 which 
directs light into single mode fiber f 2 . The interference be- 
tween this local oscillator and the scattered electric field 
is achieved by a single mode coupling device C which 
connects fibers f x and f 2 . This allows for great flexibil- 
ity when choosing the scattering angle. From the output 
of C , the electric field arising from the interference prop- 
agates through single mode fiber / 4 to a photomultiplier 
tube PUT (Thorn EMI Type 9868), which is connected 
to an electronic correlator (Malvern Series 7032 with 64 
channels). 

Between SF 2 and beam splitter BS, an attenuator A 
avoids a too high illumination of the photomultiplier tube. 
A retardation plate (A/2 plate) permits us to change the 
polarization of the light entering f 2 . The axis of the A/2 
plate is tuned in such a way that the two electric fields 



entering C have the same polarization. Thus, the inter- 
ference between the scattered electric field and the local 
oscillator measured with the PMT is optimal. 

To access the geometrical features of the scattering 
volume V, we can turn on a second laser La 2 (He-Ne, 
10 mW), which is coupled to single mode fiber / 3 via spa- 
tial filter SF 3 . Light from La 2 propagates through f 3 , then 
through f ± and goes out via SF 1 . Therefore we get a direct 
visualization of the scattering volume V as the intersection 
of the incident beam with the beam from La-,. 



4.2 Design of the Couette cell and refraction effects 




R 2 
Rs 
R 4 

Rr 



Fig. 5. Experimental Couette cell. Plexiglas (optical index n p ) 
is indicated in gray. The sample (optical index n) lies between 
Ri and R 2 and the thermostated water (optical index n w ) be- 
tween R 3 and R 4 . 



In the experimental setup described in Sec. one of 
the major constraints is the use of a thermostated Couette 
cell. Indeed, our home-made cell sketched in Fig. || is com- 
posed of four cylindrical interfaces. These interfaces sepa- 
rate the fluid under study from the water used for temper- 
ature control. This Couette cell thus acts as a cylindrical 
lens which leads to various optical complications: (i) the 
angle imposed by the operator differs from the actual 
scattering angle 8, and (ii) the scattering volume lies at 
point A instead of point A' (see Fig. |J). This leads to intri- 
cacies in the experimental determination of the positions 
of the stator and the rotor. Moreover when the Couette 
cell is moved in the x direction by a quantity Sx, the real 
displacement Sx' of the scattering volume slightly differs 
from Sx. 

All these problems due to refraction must be taken into 
account in order to get reliable velocity profiles. To explicit 
the different relationships between 9, Qi, Sx, and Sx', we 
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simply write the Sncll-Dcscartes laws for the four inter- 
faces. The geometrical parameters are as follows: R x = 
22, R 2 = 25, R 3 = 27.6, R 4 = 30.5 and R 5 = 35.6 mm. 
Hence, the aspect ratio of our Couette cell is: (R 2 — R±)/ R 2 = 
0.12. All the fluids studied in this paper have the same 
optical index n = 1.40. The optical indices of Plexiglas 
and water are n p = 1.47 and n w = 1.33 respectively. The 
imposed angle is 9i = 35°. With such parameters, Snell- 
Descartes laws lead to Sx' = f x Sx with f x ps 1.13 and 
9 = fg6i with fg ps 0.79. We checked that f x does not 
depend on the position in the gap up to a very good ap- 
proximation. However 9 depends slightly on x via fg: its 
relative variation is 69/ 9 ps 3 % across the cell gap. 

This last point is very important because the exper- 
imentally measured quantity, the Doppler shift qv, de- 
pends on the scattering angle. When scanning the width 
of the gap during an experiment, this small variation of 9 
leads to the following intrinsic incertitude on the velocity 
v: 5v/v ps 3 %. Let us notice however that this uncertainty 
strongly depends on the gap width: in a narrower gap, for 
instance e = 1 mm, Sv/v ps 1 %. 

In the next paragraph we show that the use of single 
mode fibers leads to a simplification of Eq. (|ll|) and we 
present measurements of the shape of the scattering vol- 
ume using laser La 2 . Such measurements will allow us to 
compute the theoretical correlation function. 



4.3 Numerical derivation of the heterodyne correlation 
function under shear 

4.3.1 A simplification due to single mode fibers 

A single mode fiber intrinsically collects coherent light 
only. Therefore the number of coherence areas detected 
by single mode fiber f 1 is one: the size of the photocath- 
ode area does not play any role in our case. This feature 
is co nsiste nt with the experimental results presented in 
Sec. 1.4. l| where we will stress the fact that the only rele- 
vant geometrical time is t^. According to this simplifica- 
tion, Eq. jll]) can be rewritten as: 



C(t) cx cos(q-v) 



(12) 



To proceed further in the integration of Eq. (|T^) , one has 
to know the exact shape of the scattering volume. 



4.3.2 Shape of the scattering volume 

The geometry of the scattering volume depends on the 
field diaphragm Df, the focal length of the lens used in 
spatial filter SF ± and the distance R between SF 1 and the 
scattering volume V. 

To determine th e siz e of V, we image V using laser Lai 
as explained in Sec. 4.1. With a sharp blade moving across 




the beam and perpendicular to the beam, we measured the 
transmitted intensity for different heights z of the blade 
(see Fig. |(a)). 



-500 

z (nm) 



Fig. 6. (a) Intensity collected on a photodiode by moving a 
sharp blade across the beam collecting the scattered light and 
visualized by turning on laser La2 (see Fig. ^|). (b) Intensity 
profile across the beam collecting the scattered light. The focal 
length of the microscope lens of spatial filter SFi is 2.9 mm 
(•), 5.5 mm (o), and 11 mm (□). 



From this measurement, the intensity profiles of the 
beam are easily obtained by differentiating the previous 
signals (see Fig. |^(b)). In Fig. |[ the focal length of spatial 
filter SFx was varied and the intensity profiles were mea- 
sured at the distance R = 13 cm from SF ± after the image 
of single mode fiber / t was focussed on the blade (see be- 
low). Such a method allows us to estimate the typical size 
c' of the scattering volume. We chose to define c' as the 
half- width of the intensity profiles at the relative intensity 
I/I = e~ 2 . As shown in Fig. ^(b), the minimum size is 
obtained for the focal length / = 11 mm: d ps 40 /xm. 
This is consistent with the geometrical magnification fac- 
tor 7 ps R/f ps 10 which yields d ps jd ps 40 /im, where 
d = 4 /xm is the core diameter of the fiber and R — 13 cm. 

To make sure that the image of single mode fiber f 1 
is located exactly on the blade, one has to measure the 
width of the beam at the position of the blade for different 
defocusing distances d FL between the lens and single mode 
fiber /j. In Fig. [?], we plotted the characteristic length 
c' measured for different distances d FL with the lens of 
focal length / = 11 mm. This curve shows a minimum at 
the distance d FL for which the image of the single mode 
fiber by the lens lies exactly on the blade. This yields the 
typical characteristic length c' = 40 (im reported above 
in Fig. ||(b). When the refraction effects induced by the 
Couette cell are taken into account, the length c displayed 
in Fig. I is given by c = f x c'/ sin(f e 9 t ) « 100 /im. 

Using the same method with laser La Xl we obtained 
the transverse dimension a of the scattering volume: a ps 
20 /zm. This characteristic size is consistent with the re- 
lation for Gaussian beams: a ps X f /(■nij a ) ps 20 /im where 
uJq ps 1 mm is the beam waist of the incident light and 
/ = 10 cm is the focal length of lens L. 



4.3.3 Numerical integration of C(t) 

We computed C(r) from Eq. ( p"2| ) using the experimental 
profile for the scattering volume and with a space dila- 
tion factor f x / sin# due to the Couette cell. This leads to 
the results displayed in Fig. ||. In Fig. ||(a) is shown the 
square of heterodyne correlation functions calculated from 
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Fig. 7. Width of the scattering volume c' measured from the 
intensity profiles for different "defocusing" distances d F L- The 
focal length of SFi is 11 mm. 



Eq. ([l2|) with v — for various shear rates. These func- 
tions are thus equal to |<7^( r )| 2 an d are also exactly equal 
to the normalized homodync correlation functions (see 
Eq. (0)). In Fig. ||(b), we plotted the half-time r 1/2 of the 
computed homodyne functions, defined by C(r 1/2 ) = 1/2, 
vs. various shear rates. All the qualitative features of ho- 
modyne DLS under shear are reproduced by Eq. (|l2[). 



For instance, the continuous line in Fig. 
Ti/2 oc 7 _1 , as expected from the physica 
above. 



i(b) shows that 
ideas discussed 





Shear rate (s ) 



Fig. 8. (a) Theoretical homodyne correlation functions com- 
puted from Eq. (Q) with v = for various shear rates 7 = 
1 s" 1 (V), 5 s" 1 (A), 10 s" 1 (+), 50 s" 1 (x) and 100 s" 1 (o). 
(b) Computed Ti/ 2 vs. shear rate. The solid line is the best fit 
by a power law: t 1/2 oc -y~° " 



mixture. The optical index of the suspension was mea- 
sured to be n = 1.40 ±0.01. The imposed scattering angle 
is 9i = 35°. This mixture is a Newtonian solution and, in a 
Couette cell of aspect ratio e/R 2 = 0.12, we expect to mea- 
sure velocity profiles that are very close to straight lines. 
In the following, we present experimental data showing 
that, even if in principle the homodync geometry allows 
us to estimate the local shear rate, it is not a suitable 
method to measure velocity profiles in a Couette cell be- 
cause of the strong influence of the shape of the scattering 
volume and because of the uncertainties on the exact po- 
sitions of the rotor and the stator. On the other hand, we 
show that the heterodyne geometry is a time- and space- 
resolved method to perform local velocimetry. 

Typical normalized correlation functions obtained un- 
der shear (7 = 10 s _1 ) are presented in Fig. || where the 
square root of the homodyne function is plotted in or- 
der to be compared to the heterodyne one. As expected 
from Eqs. ([!]) and (Q), C(r) 1 / 2 in homodyne mode corre- 
sponds to the envelope of the oscillating heterodyne func- 
tion. Moreover this envelope shows a strong decorrelation 
at time t « 80 /xs. 



u 




-0.5 



x (s) 

Fig. 9. Experimental homodyne (o) and heterodyne (•) corre- 
lation functions recorded on a latex suspension under shear at 
7 = 10 s _1 . The measurements were taken at xt = 14.60 mm 
i.e. 1.20 mm away from the moving inner wall. The thick line 
shows the portion of the heterodyne function that is used in the 
data analysis to determine the zero crossings (x) from which 
the local velocity is calculated (see text). 



4.4 Calibration of the experiment using a Newtonian 
fluid 

To proceed further, one has to perform calibration experi- 
ments with a well-known flow. We used a dilute suspension 
of latex spheres (<j> sa 0.1 % wt.) in a 1:1 watcr-glycerol 



4.4.1 Experimental features of homodyne DLS under shear 

Let us first consider homodyne correlation functions. In 
Fig. [lO], the half-time t 1/2 of homodyne correlation func- 
tions is plotted against the imposed shear rate 7. This 
time scales as r 1/2 oc j^ 1 . Since the characteristic geo- 
metrical times scale as 77 oc 7 _1 and T4 oc 7 _1 , this is 
the scaling expected from theory. Moreover the measured 
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Fig. 10. Half-time Ti/ 2 of the homodyne correlation function 
measured in the middle of the gap for (Ifl — 200 /im vs. shear 
rate. For each data point, the correlation function was accumu- 
lated over 1 min. The solid line is a power-law fit with exponent 
-0.99. 



Fig. 11. Half-time Ti/2 of the homodyne correlation function 
measured in the middle of the gap at 7 = 10 s _1 vs. "de- 
focusing" distance cLfl- For each data point, two correlation 
functions were accumulated over 1 min and the two values of 
T1/2 were averaged. The error bars account for the difference 
between those two estimates. 



Ti/2 (7) are in good agreement with those computed from 
Eq. (||) and displayed in Fig. |(b). 

In Fig. |l|, we plotted r 1/2 for different defocusing dis- 
tances d FL at a given shear rate 7 = 10 s _1 . As explained 
above, changing the distance d FL between the lens and the 
single mode fiber increases the width of the beam at a fixed 
distance R. Therefore one expects an increase of the size of 
the scattering volume V. As shown in Fig. Rl, this leads to 
a decrease of the characteristic time r 1/2 . This behaviour 
can be expected from theory if one considers oc (erf) . 
From these results, it seems that the relevant time which 
governs the correlation function C(t) is r^. It is in agree- 



ment with the features discussed in Sec. 4.3, where we 
emphasized that should not play any significant role 
thanks to the use of single mode fibers. 

Homodyne DLS allows quantitative measurements of 
the local shear rate only if the shape of the scattering 
volume V is precisely known. In our experiment, the scat- 
tering volume is obviously not well defined mainly because 
of the refraction effects induced by the Couette cell. For 
instance a small motion the Couette cell can change the 
shape of the scattering volume. 

In Fig. |l|, we plotted the half-time r 1/2 of homodyne 
correlation functions for various positions Xt of the me- 
chanical table and for two different defocusing distances 
d FL . For d FL = 0, the values of r 1/2 show a dispersion of 
about 20 % and for d FL = 200 /im, the dispersion is about 
50 %. In both cases the half-times significantly grow when 
approaching the two walls. This effect, which accounts for 
most of the observed dispersion, is due to the fact that 
the size of scattering volume decreases when V intersects 
one of the cell walls. This increase of r 1/2 is particularly 
noticeable when d FL = 200 because in this case the 
scattering volume is about twice as large as for d FL = 
(see Fig.0). 



These intricacies induced by the Couette geometry do 
not allow for simple quantitative measurements of the lo- 
cal shear rate. Moreover with such measurements, we are 
not able to determine the exact positions of the rotor and 
the stator of the Couette cell. These are the reasons why 
we chose heterodyne DLS for local measurements under 
shear. 
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Fig. 12. Half-time T1/2 of the homodyne correlation function 
vs. position in the gap for a well-focussed collection optics 
dFL = (o) and for a defocussed case dFL — 200 /im (•). 
For each data point, the correlation function was accumulated 
over 1 min. The shear rate was set to 7 = 10 s _1 . 
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4.4.2 Experimental features of heterodyne DLS under shear with R 1 and R 2 denoting the radii of the rotor and the 



In order to extract the Doppler shift q • v from the ex- 
perimental heterodyne correlation function, we interpo- 
late the experimental curves as shown in Fig. ^. A series 
of cancellation times Tk is extracted from the interpola- 
tion. The Tfc's arc marked by x symbols in Fig. ^. Since 
C(t) oc F(t) cos(q-vr). one expects = (/c7r+ : |)(q-v) _:L 
provided F(t) does not vanish. It is then straightforward 
to estimate the scalar product qv from a linear fit of 
vs. k. 

To obtain velocity profiles, one has to move the scat- 
tering volume across the gap of the Couette cell. This is 
done by moving the rhcometer in the x direction with the 
help of the mechanical table (cf . Fig. ||) . The Doppler pe- 
riod T = (qv) -1 is plotted for different positions Xt of the 
table and for different imposed shear rates in Fig. |l3|(a). 
To estimate the velocity corresponding to a given Doppler 
period, one has to k now precisely the scattering vector q. 
As explained in Sec. 4.2 refraction leads to important op- 



tical corrections. Indeed the real scattering angle 9 differs 
from the imposed angle 0i since « feOi, where fg w 0.79. 
Such a correction allows us to convert the measured period 
T into a velocity v as shown in Fig. |l3|(b). Indeed, assum- 
ing that v = ve y , i.e. that the flow is purely orthoradial 
in the (x, y) plane, q-v = qv cos(9/2) = 2ivnsm6v/X. 
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Fig. 13. Calibration of the velocity measurements using a 
Newtonian suspension of latex spheres in a water-glycerol mix- 
ture submitted to different shear rates 7 = 2 s _1 (•), 5 s _1 (o), 
10 s _1 (x), and 20 s" 1 (□). (a) Raw estimates of the period 
T of the oscillations in the heterodyne correlation function as 
a function of the position of the mechanical table x t . The ori- 
gin of xt is arbitrary and is defined when the actuators are 
first initialized, (b) V elocity estimates using the procedure de- 
scribed in Sec. 4.4.2. Each data set is fitted by a straight line, 



stator of the Couette cell. Fa, is a conversion factor inher- 
ent to the rheometer that account for the Couette geom- 
etry. The velocity of the mobile inner wall is thus given 
by: v = R x f2 = R^/Fa,. This permits the determination 
of the position of the moving inner wall x m = 13.49 mm 
(left dashed line in Fig. |l3|(b)). Such a determination leads 
to an effective gap width e' = Xf — x m m 2.65 mm, 
which differs from the real gap e = 3 mm. Actually as 
shown in Sec. 4.2, the actual displacement of the scatter- 



ing volume inside the cell Sx' is related to the displace- 
ment of the mechanical table 5x, by the following rela- 
tion: Sx' = f x Sx, where f x ~ 1.13 is deduced from the 
Sncll-Dcscartcs law for refraction. Such a conversion fac- 
tor is in excellent agreement with the previous result since 
e/e' 1.13. Let us notice that for all the imposed shear 
rates, we were able to measure velocities for a few points 
that are outside the gap at xt < x m . (see Fig. |l3|(b)). This 
is due to the finite extent of the scattering volume: even 
for xt < x m , some small intensity may be detected that 
is scattered from the intersection of the scattering volume 
and the gap. 

In the velocity profiles displayed in Fig. the two 
optical conversion factors fg and f x have been taken into 
account and x denotes the real position inside the gap. 
In all cases, the uncertainty of the experimental measure- 
ments is of the order of the marker size. Note that the 
effect of non-homogeneous shear rate due to the curved 
geometry of the Couette cell is negligible for such a New- 
tonian fluid and docs not influence the calibration. 




which yields an estimate for the position of the fixed outer 
wall: Xf — 16.15 ± 0.02 mm (dotted line). The dotted line 
at x m — 13.49 mm shows the estimated position of the mov- 
ing inner wall. The solid lines are the linear velocity profiles 

v = (Rij/F,,) (xt - xt)/(x f - x m ) with 7 = 2, 5, 10, and Fi §- 14 - Velocity profiles in a sheared Newtonian suspension 

20s -l. for 7 = 2 s- 1 (.), 5 s" 1 (o), 10 s" 1 (x), and 20 s" 1 (□). The 

abscissa x denotes the real position within the cell gap. The 
origin of x is taken at the inner wall. The solid lines are the 
Linear fits of the measured velocity profiles allow us linear velocity profiles v = R^l-x/ej/F^ with 7 = 2, 5, 10, 

to estimate the position of the fixed outer wall Xf = and 20 s" 1 . The conversion factors f g — 0.79 and f x = 1.13 

16.15 mm with a precision of about 20 ^m. The velocity v 

of the rotor is given by the rheometer: for a given angular 

velocity Q of the rotor, the shear rate 7 indicated by the 

rheometer, 7 = F^fl where Fa, = (R* + Rl)/{R 2 2 - R%) 



are used to account for refraction effects 
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4.4.3 Spatial resolution 

One of the most relevant features of a technique measur- 
ing velocity profiles is its spatial resolution. As previously 
outlined, a good resolution would be around 100 fim. In 
our case, spatial resolution is limited by the size of the 
scattering volume V. Indeed, velocity measurements using 
heterodyne DLS are averaged over V. The typical dimen- 
sions of V in air are a w 20 /xm and c! f» 40 /im mea- 



to follow dynamics like those observed in Rcfs. [|11|,|12| that 
involve time scales in the range 1-10 3 s. 



sured in Sec. 4.3.2. Since scattering occurs at an angle 
6 and since the Couette cell acts lens, the relevant 
size is c w f x c' / sm(fg9i), which leads to a typical size 
c ps 100 fim. However, the spatial resolution of our setup 
strongly depends on the scattering angle. The latest mea- 
surements performed in our group on banded flows in mi- 
cellar solution at a larger angle (6*; r* 50°) reveal that our 
setup can resolve a discontinuity of the shear rate with an 
accuracy of about 50 /im p9| . 



4.4.4 Temporal resolution 

As recalled in the Introduction, some complex fluids dis- 
play non-stationnary flows [{ll[ , |l2| . Therefore it is inter- 
esting to perform time-resolved measurements of the local 
velocity. To obtain the Dopplcr frequency with a high tem- 
poral resolution one may use a spectrum analyzer rather 
than a correlator. Such a technique allows up to a thou- 
sand measurements per second. For instance, this has been 
implemented long ago by Gollub et al. to measure the tem- 
poral evolution of the local velocity in Newtonian fluids at 
the onset of the Taylor-Couette instability |2(|. However 
the use of a correlator yields more information than a spec- 
trum analyzer. In homodyne mode qualitative information 
on the local shea r rate can be obtained, as discussed in 
Sees. S and 4.4.1. Moreover, the use of a correlator once 



the shear flow has been stopped can give important infor- 
mation on the coupling between the structural dynamics 
of a material and an applied shear || . 

The temporal resolution of our setup strongly depends 
on the power of the laser beam and on the quantity of 
light scattered by the sample. Actually these two features 
control N(t), the number of photons detected by the pho- 
tomultiplier tube per second. To estimate the temporal 
resolution of our setup, we measured the homodyne and 
heterodyne correlation functions for various accumulation 
times t a at a given N(t). Then we computed the corre- 
sponding half-times of the homodyne correlation functions 
and the Doppler periods in the heterodyne mode. Fig- 
ure 15 displays such measurements obtained on a latex 
suspension in the middle of the gap for 7 = 10 s _1 and 
with N(t) w 40 4 events per second in homodyne mode. 
The error bars represent the standard deviations of differ- 
ent measurements obtained at a given accumulation time 
t a . It clearly appears that good statistical estimates of r 1/2 
and T = (q-v) -1 are obtained after an accumulation time 
of the order of 5 seconds. Actually in better experimental 
conditions, for instance when N(t) w 3.10 4 , the accumu- 
lation time required to define those characteristic times 
decreases to 1 second. Therefore our setup is well suited 



1 m\ 



(b) 



lo s,„(t) 



Fig. 15. Statistical convergence of the estimates of the charac- 
teristic times, (a) Half-time T1/2 of the homodyne correlation 
function vs. accumulation time t a . (b) Period T of the oscilla- 
tions in the heterodyne correlation function vs. accumulation 
time t a - Each data point represents the average over the the 
characteristic times estimated from five consecutive correlation 
functions accumulated over t a . The error bars account for the 
standard deviation of those five measurements. The shear rate 
was set to 7 = 10 s _1 . 



5 Discussion and conclusion 

In this article, we presented an experimental setup based 
on the use of single mode fibers that allows one to mea- 
sure velocity profiles in Couette flows. We have shown 
that measurements of the local shear rate using homo- 
dyne DLS are difficult to implement in the Couette ge- 
ometry. By measuring directly the shape of the scattering 
volume, we were able to compare experimental and theo- 
retical homodyne correlation functions, and show that the 
decorrelation times t 1/2 are linked to the size of the scat- 
tering volume c and to the local shear rate 7 according to 
Ti/2 oc 1/(07) as expected from theory. 

Heterodyne detection yields the Doppler shift q • v as- 
sociated to the average velocity in the scattering volume. 
Provided refraction effects are taken into account and a 
careful calibration is performed, our setup gives access to 
velocity profiles in Couette flows of complex fluids down 
to a spatial resolution of 50 /im. As far as temporal resolu- 
tion is concerned, the major limitation is the time required 
to move the rheometer along the direction of shear. In the 
best experimental conditions and with about 30 points in 
the gap of the Couette cell, we managed to measure veloc- 
ity profiles in about one minute. Thus, our setup cannot 
resolve flow dynamics over the whole gap on time scales 
shorter than one minute. Other non intrusive techniques 
like ultrasonic speckle correlation recently developed in 
our group may allow us to measure up to 10 2 velocity pro- 
files per second and may yield rich information on faster 
dynamics. 

In order to show its relevance in the field of complex 
fluid flows, the DLS technique described in this paper was 
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Fig. 16. Velocity profiles for 7 = 2 s _1 in two different emul- 
sions, (a) Dilute emulsion <f> = 20 %. (b) Concentrated emulsion 
4> = 75 %. The dashed lines represent the velocity profiles ex- 
pected for a Newtonian fluid. The dotted lines correspond to 
a Newtonian fluid when wall slip has been taken into account 



applied to sheared emulsions. Figure [l6| shows the veloc- 
ity profiles measured in two different oil-in-water emul- 
sions at the same shear rate. The Couette cell is that used 
throughout this study (e = 3 mm and i? x = 22 mm). Both 
emulsions are composed of silicone oil droplets (diameter 
« 2 /mi) in a water-glycerol mixture with 1 % TTAB sur- 
factant. Figure [l6|(a) corresponds to a dilute emulsion of 
volume fraction cf> = 20 % whereas Fig. |l6|(b) was obtained 
with a concentrated system of volume fraction <fi = 75 %. 
For 7 = 2 s" 1 . rhcological data recorded by the rheometer 
show that the global stresses imposed on the samples are 
a = 0.07 Pa and a = 290 Pa respectively. The optical 
index n of the aqueous phase is matched to that of the sil- 
icone oil to avoid multiple scattering. Moreover, n = 1.40 
as in the latex suspension used for the ca libra tio n, so t hat 
the values of f x and fg discussed in Sec. 4.2 and 4.4.2 re- 



main valid in the case of those emulsions. In both cases, 
slippage clearly shows on the velocity profiles when com- 
pared to that of a Newtonian fluid. Wall slip is much larger 
for the concentrated system. Moreover, the velocity pro- 
file is almost perfectly linear in the dilute emulsion (see 
fig. |l^(a)) whereas a small but systematic curvature can 
be noticed in the concentrated emulsion (see fig. p~6|(b) ) . 
This curvature is the signature of the non-Newtonian be- 
havior of the 75 % emulsion. The complete analysis of the 
slip velocities and of the bulk flow behavior inferred from 
such local measurements will be discussed at length in a 
related publication pc|| . 

Finally, the main innovation of the present work is 
to provide both local velocity measurements and global 
rheological data simultaneously. DLS coupled to standard 
rheology opens the door to local rheology. Indeed, when 
the flow is stationnary, the local stress in a Couette cell 
is inversely proportional to the square of the position in 
the gap. On the other hand, the local shear rate may be 
accessed by differentiating the velocity profiles. Thus, we 
could determine the relationship between the local shear 
rate and the local stress and define a local flow behaviour. 
The comparison between such a local rheology and the 
global flow curve obtained with a standard rheometer sho- 
uld lead to essential knowledge that ma y p rove crucial to 
understand complex fluid flows 0,^,LM[jO| . 
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